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Abstract:  This  paper  investigates  the  stability  properties  of  second  order 

systems,  x  =  f(x),  where  f(x)  contains  either  quadratic  terms  -system  (1)-  or 
linear  and  quadratic  terms  -system  (2)-in  x.  The  principal  contributions 
are  summarized  in  two  theorems  which  give  necessary  and  sufficient  conditions 
for  stability  and  asymptotic  stability  in  the  large  of  systems  (1)  and  (2) 
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The  Stability  of  Second  Order 
Quadratic  Differential  Equations 

Daniel  E.  Koditschek  and  Kumpati  S.  Narendra 


1 .  Introduction: 

The  diverse  behavior  of  nonlinear  systems  reflects  the  obvious  fact  that  no  one 

theory  can  address  all  phenomena  of  interest.  Powerful  techniques  have  evolved  over 

the  years  in  the  context  of  specific  questions  regarding  particular  classes  of  such 

systems.  This  paper  suggests  a  useful  approach  to  the  investigation  of  the  stability 

characteristics  of  a  class  of  second  order  differential  equations: 

T 

x  Gx 

T  =  B(x)  (1) 

x  '  ;  x  Hx  1 

L  2J  L 

and 

x  =  Ax  +  B (x)  (2) 

More  than  a  convenient  choice  of  systems,  quadratic  differential  equations  have 

a  traditional  and  increasing  importance  in  stability  literature.  Given  a  general 
autonomous  differential  equation  x  =  f(x),  f(x)  can  be  expanded  in  a  Taylor  Series 
around  the  origin  if  it  satisfies  certain  regularity  conditions.  If  the  equilibrium 
state  of  the  linear  approximation  is  asymptotically  stable  (unstable)  it  is  well- 
known  that  the  nonlinear  system  is  also  asymptotical ly  stable  (unstable).  However, 
in  the  critical  case,  when  the  linear  approximation  is  merely  stable,  interest  shifts 
to  an  equation  of  the  form  (2)  in  the  stability  analysis  of  the  nonlinear  system. 
Differential  equations  of  this  type  are  known  to  arise  in  adaptive  control  where 
control  parameters  become  state  variables  of  a  quadratic  system.  More  recently,  the 
special  class  of  bilinear  systems 


x  =  Ax  +  uDx  +  bu 


(3) 
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has  received  a  great  deal  of  attention  in  the  control  literature  and  the  principal 
results  of  this  theory  are  influencing  the  direction  of  research  on  general  nonlinear 
systems.  When  the  control  u(t)  in  equation  (3)  is  a  linear  function  of  the  state 
variables,  equation  (3)  becomes  a  special  case  of  (2)  and  the  stability  properties 
of  such  systems  are  bound  to  be  of  interest  to  control  theorists. 

This  paper  presents  necessary  and  sufficient  conditions  for  stability  in  the 
large  of  system  (1)  and  asymptotic  stability  in  the  large  of  system  (2).  A  brief 
treatment  of  some  interesting  properties  of  general  homogeneous  differential 
equations  in  section  2  provides  the  setting  for  the  results  on  pure  quadratic 
systems  developed  in  section  3.  A  simple  statement  of  conditions  for  the  stability 
of  system  (1)  is  given  by  Theorem  1  at  the  end  of  that  section.  While  the  quali¬ 
tative  behavior  of  (1)  can  be  completely  characterized  using  a  few  simple  prototypes 
[  1]  the  addition  of  the  linear  part  in  (2)  results  in  much  more  varied  behavior: 
a  simple  classification  will  not  be  complete;  a  complete  classification  will  not 
be  simple.  Consequently,  we  limit  our  interest  to  one  important  aspect  of  the 
stability  behavior  of  (2)  -  asymptotic  stability  in  the  large  (a.s.i.).  Yet  even 
within  this  restricted  scope,  only  necessary  or  sufficient  conditions  for  a.s.£. 
are  amenable  to  simple  presentation.  A  complete  account  is  seen  in  section  4  to 
involve  a  variety  of  special  cases.  The  resulting  statement  of  necessary  and 
sufficient  conditions  for  a.s.H.  is  given  by  Theorem  2  at  the  end  of  that  section. 

2 .  Homogeneous  Systems  of  Even  Degree: 

n 

Consider  a  dynamical  system  in  R  : 

x  =  h(x)  (4) 

where  h(x)  is  a  homogeneous  function  such  that  h(8x)  =  bVx),  6  e  R  k  >  1.  If 
k  is  even  the  direction  of  the  field  is  a  constant  along  any  straight  line  through 


the  origin.  A  well-known  but  useful  fact  concerning  solutions  of  (4)  is  given  by 
the  following  lemma. 

Lemma  2.1:  Let  p(t;x^)  be  a  solution  of  (4)  given  ^ 
initial  condition  p(0;Xg)  =  x^.  Then  for  all  8  e  IR  S' 


k-1 

=  Bp  (8  t;xQ) 

Proof •  Let  t  =  g  s  and  define  u(s)  =  pp [ t ( s) ; ] . 
Then  ^u(s)  =  ^p(t;x0)g=8k^p(t;x0) 


=  h(u) 


Figure  1.  Reflection  Property  of  Even 
Homogeneous  System 


Hence  u(s)  satisfies  (4)  with  initial  condition  u(0)  =  8xq  as  does  p(t;BxQ)  which 
implies  u(s)  =  p(t;8xQ).  Q 

If  k  is  even  then  this  relationship  implies  p(t;-x^)  =  -p(-t;xQ).  In  other 
words  any  trajectory  through  xQ  for  t  ^  0  has  an  associated  trajectory  through 
~XQ  for  t  j:  0  which  is  its  reflection.  This  simple  fact,  depicted  in  Figure  1, 
leads  to  the  following  corollary. 

Corollary  2.1:  If  k  is  even  and  the  origin  is  stable  then  for  any  x^  /  0  the 
complete  trajectory  y(Xq)  =  {p(t;x^)|t  z  01}  is  a  positive  distance  from  the  origin. 
Proof :  Let  ||xj|  >  e  >  0  but  let  the  positive  trajectory  not  be  bounded  away  from 
the  origin.  Then  there  exists  a  sequence  { t  }  ,t  >0  such  that  ||  p ( t  ;x_)(l  <  — 

n  n  II  n  0  II  n 

Defining  x  =  p(t  ,x  )  we  have  p(t  ;-x  )  =  -x„  from  Lemma  2.1.  Hence  for  no  n  >  0 
n  n  0  K  n  n  0 

does  ||xn|  <  —  imply  |jp(t;-x^)|  <  e  and  the  origin  is  unstable.  A  similar  argument 
applies  to  the  negative  trajectory.  □ 


In  general.  Corollary  2.1  indicates  that  an  even  degree  homogeneous  dynamical 

n  2 

system  in  R  can  never  be  asymptotically  stable.  On  the  plane,  R  ,  we  may  say 
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even  more.  Equation  (1)  is  a  specific  example  of  (4)  in  R  for  which  k  =  2. 

2 

The  remainder  of  this  paper  will  be  devoted  to  the  study  of  (1)  and  (2)  in  R  . 

Lemma  2.2:  No  solution  of  system  (1)  other  than  an  equilibrium  point  can  be  a 
closed  path. 

Proof :  This  follows  directly  from  the  fact  that  any  nontrivial  closed  curve  will 

intersect  some  line  through  the  origin  at  least  twice.  A  trajectory  on  this  curve 
would  imply  that  the  field  changes  direction  along  that  line  which  violates  the 
even  homogeneous  property  of  (1).  □ 

Corollary  2.2:  If  the  origin  of  system  (1)  is  stable  then  the  field  must  vanish 
along  at  least  an  entire  line  through  the  origin. 

Proof :  It  suffices  to  show  that  an  equilibrium  state  t  0  exists.  By  homogeneity 

the  conclusion  follows. 

Let  V  ?  0  be  a  trajectory  of  (1)  contained  in  some  compact  neighborhood  of  0. 
From  Corollary  2.1  0  t  y  where  y  is  the  closure  of  y .  From  Lemma  2.2  Y  is  not  a 
non-trivial  closed  path.  Hence  by  Poincare-Bendixson  Theorem  the  limit  set  of  y  is 
a  singular  point  xQ  of  B.  Since  xQ  t  0  we  have  B(x)  =  0  for  all  x  =  ax^,  a  e  R 
by  homogeneity.  □ 

The  next  section  will  investigate  the  existence  of  lines  along  which  the 

field  in  (1)  vanishes  as  required  by  the  previous  corollary,  leading  to  a  new 

parametrization  of  stable  quadratic  systems.  This  new  parametr ization  allows 

2x2 

the  characterization  of  stability  behavior  in  terms  of  a  matrix  in  R  ,  a  result 
whose  consequences  pervade  the  remainder  of  this  paper.  Indeed,  concerning  the  final 
section  of  the  paper,  the  availability  of  this  linear  form  can  safely  be  said  to  make 
a  seemingly  hopeless  task  complicated  but  possible. 
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3.  Stability  of  Second  Order  Quadratic  Systems: 

In  this  section  we  shall  exclusively  consider  the  particular  class  of  second 
degree  second  order  systems  described  by  (1) 


T 

x  Gx 


=  B(x) 


(1) 


We  assume  that  at  least  one  of  either  G  or  H  e  R  is  non-zero  and,  without  loss 
of  generality,  that  both  are  symmetric.  From  Corollary  2.2  it  follows  that  the 
locus  of  the  set  of  critical  points  of  B(x)  is  crucial  to  the  stability  properties 
of  (1).  This  is  completely  determined  by  G  and  H,  and  using  familiar  properties  of 
symmetric  matrices  we  may  classify  the  types  of  equilibrium  states  evinced  by  (1) 
according  to  whether  the  field  vanishes. 

(i)  only  at  the  origin 

(ii)  along  a  straight  line  through  the  origin 

(iii)  along  two  straight  lines  through  the  origin. 

Table  I  provides  examples  of  each  type  of  quadratic  field  based  upon  the  familiar 
rank  and  signature  properties  of  G  and  H.  The  reader  may  note  that  as  an  immediate 
consequence  of  Corollary  2.2,  systems  of  type  ( i)  cannot  be  stable  -  in  the  terms 
of  Table  I  only  systems  like  those  given  by  examples  2b,  2c,  3b,  or  4b  may  be  stable. 
The  remainder  of  this  section  is  devoted  to  refining  this  necessary  condition, 
a .  Notation  and  Definitions: 

Since  indefinite  and  semi-definite  matrices  will  arise  in  all  subsequent  dis¬ 
cussions  it  is  worth  establishing  the  following  notational  conventions  concerning 
their  algebraic  and  geometric  properties.  Let  Ag  denote  the  symmetric  part  of  A. 

The  set  S [ A ]  =  {M|a  =  M^}  denotes  the  symmetric  equivalence  class  of  A.  Clearly, 

T  T  2  T 

x  Mx  =  x  Ax  V  x  e  R  if  and  only  if  M  e  SlA].  If  A  =  ab  is  singular  and  P  e  S[A] 

T 


then  x  Px  =  0  iff  x  is  orthogonal  to  either  a  or  b. 


table  I 


Example 


G  or  H  definite 

Type  (i) 

a)  G  and  H  Indefinite 

Type  (i) 

and  of  full  rank. 

b)  G  and  H  indefinite 

Type  (ii) 

and  of  full  rank. 

c)  G  and  H  indefinite 

Type  (iii) 

and  of  full  rank.  j 

i 

! 

G  singular  H  indefinite  i 

and  of  full  rank. 

a)  Gx  =  0  xTHx  /  0 

j 

1 

i  Type  (1) 

T  1 

b)  Gx  =  0  x  Hx  =  0 

Type  (ii) 

G  and  H  singular 

a)  G  ^  aH 

Type  (i) 

b)  G  =  aH 

Type  (ii) 

i 

0  1/2 ~ 
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The  class  S[0]  contains  a  matrix  J  with  the  property  J  =  -I .  J  Is  the 


skew-symmetric  matrix 


0  -1 


1  0 


0  1 


-1  0 


and  maps  every  vector  In  1R  into 


its  orthogonal  complement  in  the  counterclockwise  (clockwise)  sense  -  i.e.  x  =  Jx. 

T  2  T 

From  these  definitions  it  follows  that  (i)  xx  =0  V  x  £  R  ;  (ii)  x  y  =0  iff 

C  „  ,  "1  i  J  i 


x  =  ay  a  e  r;  (iii)  For  J  =  x  y 

1  oj  1 

denotes  the  determinant  of  the  array  (y,x). 


0  -1 


T  T 

x  =  x  Jy  =  |y,x|  where  the  last  symbol 


Finally,  the  following  conventions  will  facilitate  further  discussion  involving 

2 

some  natural  geometric  objects  defined  by  vectors  and  matrices  in  R  .  We  denote 

the  subspace  defined  by  x  £  R2  as  | (x)  =  (y  |y  =  ax,  a  e  R),  its  orthogonal  complement 

as  i(x)  =  (y  |  yTjc  =  0}  and  its  positive  affine  half  plane  as  H^(x)  =  {y|xTy  >  q} 

for  b  £  R.  If  P  is  a  symmetric  indefinite  or  semi-definite  matrix  then  its  positive 

A  i  T  2 

cone  is  defined  as  C  (P)  =  (x|x  Px  >  y}  for  y  e  R.  If  [)  c  R  is  disconnected  [4, 
p.  152]  then  there  exists  a  pair  of  open  sets  0^,0^  such  that  0  =  0^u02  and  0  =  $ 

called  a  separation  of  0.  The  set  0^  is  said  to  separate  0  and  separates  improperly 
when  0^  =  j>. 

b.  Quadratic  Differential  Equations  with  Non-Isolated  Equilibria: 


These  considerations  allow  an  immediate  re-parametrization  of  system  (1)  with 
non-unique  equilibria. 

2 

Lemma  3.1:  System  (1)  is  type  (ii)  or  (iii)  if  and  only  if  there  is  some  c  e  R 

2x2  T 

and  D  £  R  such  that  B(x)  =  c  xDx. 


Proof :  B(x)  is  type  (ii)  or  (iii)  if  and  only  if  for  some  x^  ¥  0 


B(x1)  = 


xiGxi 

x["*l 


T  T 

Now  xjGx^  =  0  if  and  only  if  G  e  Sfd^c  ]  where  c  c  i(x^)  and  is  an  arbitrary  point 

2  T  T 

in  R  •  Similarly  H  e  Sfd„c  ].  Then,  defining  .  d..  we  have 

D  =  -- — 


T  T 

x  c  d 

T  T 

x  c  d^X 


we  have 


=  c  xDx. 
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The  importance  of  this  lemma  is  that  it  relates  individual 
solutions  of  a  nonlinear  time  invariant  system  to  solutions  of  a  linear  time  varying 
system 

x  =  Y ( t) Dx 


which,  in  turn,  may  be  written  as 


^  =  Dx(s) 
ds 


=  Y(t) 


Thus,  if  p(t;x^)  is  a  solution  to  system  (1)  of  type  (ii)  or  (iii)  then 


,  .  s(t)D 

Pit  JXq)  =  e  xQ 


s(t)  =  J  c  p(i;Xg) dx 

0 


and  that  trajectory  remains  on  the  integral  curve  of  system  (5)  specified  by  x^ 

with  position  parametrized  by  s(t).  The  following  classification  of  system  (5) 

provides  the  framework  for  the  complete  stability  analysis  of  (1)  and  the  a.s.i. 

* 

analysis  of  (2).  We  will  say  D  is  focal  if  it  has  two  complex  conjugate  eigen- 

*  * 
values,  nodal  if  it  has  two  distinct  real  eigenvalues,  and  critical  if  it  has  a 

unique  real  eigenvalue  of  multiplicity  two.  We  further  distinguish  among  the  class 

of  critical  matrices  those  which  admit  only  one  eigenvector  as  opposed  to  an  in- 

2 

finite  number  of  eigenvectors.  If  D  is  critical  with  the  unique  eigenvector  x  e  R 
then  we  will  say  D  is  x-critical.  To  conclude  this  section  we  show  that  only  focal 
matrices  admit  a  stable  origin  of  system  (1). 
c.  Necessary  and  Sufficient  Conditions  for  Stability: 


To  exclude  all  critical  and  nodal  D  matrices  from  further  consideration  we  apply 

Chetaev's  Instability  Theorem  [  1  1  as  follows. 

T 

Lemma  3.2:  If  B(x)  =  c  xDx  and  D  is  nodal  or  critical  then  system  (1)  is  unstable. 

T 

Proof :  Since  B(x)=  (-c  x) (-Dx)  we  assume  with  no  loss  of  generality  that  D  has  an 

eigenvalue  in  R  .  Hence  by  Fact  A. 4  in  Appendix  A,  there  exists  an  indefinite 

* 

Note  the  distinction  between  these  classes  and  the  traditional  designation  of 
linear  systems  according  to  whether  the  origin  is  a  node,  saddle,  focus,  or  center. 
For  instance,  a  focal  matrix  defines  a  linear  system  whose  origin  Is  either  a  focus 
or  center,  while  a  nodal  matrix  yields  a  node  or  saddle. 
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symmetric  matrix  P  such  that  c.  t  C  (P)  and  C  (P)  c  ^„([PD]  ).  From  Fact  A. 3  we 

r  U  U  -  U  s 

have  il  =  ^(OnC^P)  separates  C  (P)  as  depicted  in  Figure  2.  Thus  if  we  define 

V(x)  =  1/2  x  Px  then  v  =  c^x  x^PDx,  and  V  [  R  ]  =  C_(P)  with  A  =  C _ ( [ PD ]  )nfL(c) 

0  0  s  U 

• -1  +  -1  + 

c  V  [  R  ].  Since  S2  separates  V  l  R  1  we  have  v  (x)  =  0  for  x  on  the  boundary  of 

.-1  + 

Pi  and  V(x)  >  0  for  x  e  Si.  Moreover,  Si  <  A  '  V  [  R  ]  so  we  have  a  region  where 

V(x)  >  0  embedded  in  a  region  where  v(x)  >  0  on  whose  boundary  V(x)  =  0  which 


If  D  is  focal  then  the  integral  curves  of  (5)  intersect  every  direction  on  the 

plane  at  a  finite  point  at  least  once.  In  this  case  all  trajectories  of  (1)  must 

tend  eventually  to  i(c)  which  must  therefore  contain  a  globally  attractive  invariant 

set.  A  simulation  of  such  a  system  is  given  in  Figure  3.  That  this  situation 

corresponds  to  the  stability  of  the  origin  is  proven  as  follows. 

T 

Lemma  3.3:  If  B(x)  =  c  xDx  and  D  is  focal  then  (1)  is  stable. 

Proof :  If  D  is  focal  then  o(D)  =  {oH-iS,a-i8}  a, 8  e  R,  8  ^  0.  Hence,  for  t  =  8/2n, 

Dt 

the  natural  frequency  of  linear  system  (6),  we  know  e  Xg  e  | (x^) .  From  section  3.b 
we  have 


.  Ds(t) 
p(t;xQ)  =  e  xQ 


ft 


;xn)di  t  0  for  x  i  i(c) 


s(t)  -  j  C  pf^;Xg^-.  r  ~  ~ .  — Q 


X  rs 

Thus,  if  for  some  t  e  R  s(t  )  =  t^  then  p(t  ;Xg)  e  | (Xg)  violating  the  homogeneity 


argument  in  section  2  which  constrains  all  solutions  of  (1)  to  lie  strictly  within 


F igure 
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one  half-plane  of  R  •  This  implies  s(t)  t  (-t^,t  )  V  t  e  R. 

A 

y  =  sup 


I  we  have 


|p(t;x0)||  <  )e 


Ds(t) 


xJ  <  Y  lx. 


Then  if 


and  the  system  is  stable.  !.] 

On  the  basis  of  Lemmas  3.1-3  we  may  state  the  principal  result  of  this  section. 

T 

Theorem  1 :  System  (1)  is  stable  in  the  large  if  and  only  if  B(x)  =  c  xDx  and  D 
is  focal. 

T 

Proof :  By  Lemma  3.1  system  (1)  is  not  stable  unless  B(x)  =  c  xDx.  By  Lemmas  3.2 
T 

and  3.3,  x  =  c  xDx  is  stable  if  and  only  if  D  is  focal.  By  Lemma  2.1  if  (1)  is 
stable  it  is  stable  in  the  large.  H 
d .  Summary: 

The  results  of  this  section  indicate  that  most  second  order  quadratic  differ¬ 
ential  equations  are  unstable:  Only  those  systems  whose  solutions  lie  on  the  inte¬ 
gral  curves  of  a  linear  system  whose  equilibrium  state  is  a  center  or  a  focus  may 
be  stable.  An  exhaustive  account  of  the  qualitative  behavior  of  the  far  more 
pervasive  unstable  examples  of  system  (1)  is  given  in  [  1  ].  Since  this  paper  is 
concerned  solely  with  questions  of  stability,  we  are  content  here  to  ignore  those 
results  and  immediately  apply  the  techniques  used  in  the  proof  of  Theorem  1  to  the 
problem  of  a.s .1.  for  system  (2). 

A .  Asymptotic  Stability  in  the  Large  of  System  (2):  x  =  Ax  +  B(x) 

As  mentioned  in  the  introduction,  the  addition  of  the  linear  part  in  equation 

(2) 

x  =  Ax  +  B(x)  (2) 


results  in  such  diversity  of  behavior  that  we  confine  our  attention  exclusively  to 
one  question  -  the  a.s.£.  of  equation  (2).  In  spite  of  this  limited  range,  the 


,X-V  '  '  "  - 
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complexity  of  the  system  is  such  that  several  special  cases  have  to  be  considered 
before  a  complete  solution  can  be  presented.  For  ease  of  exposition  necessary 
conditions  are  derived  in  section  4.b  and  conservative  sufficient  conditions  are 
stated  in  section  4.c.  Consideration  of  special  cases  in  4.d  to  which  the  results 
of  sections  4.b  and  4.c  do  not  apply  permit  the  proof  of  Theorem  2  which  contains 
necessary  and  sufficient  conditions  for  the  a.s.Z..  of  system  (2). 
a.  Preliminary  Discussion: 

The  existence  of  a  positive-invariant  set  disconnected  from  the  origin  is  a 
sufficient  condition  for  the  equilibrium  state  of  system  (2)  not  to  be  globally 
asymptotically  stable.  A  systematic  study  of  conditions  for  the  existence  of  such 
sets,  using  Propositions  4.1  and  4.2,  leads  to  the  necessary  conditions  of  section 
4.b  as  well  as  allowing  the  classification  of  special  cases  in  section  4.d. 
Proposition  4.1  concerns  the  existence  of  an  equilibrium  state  removed  from  the 
origin  (hereafter  referred  to  as  an  off-origin  equilibrium),  while  Proposition  4.2 
introduces  a  device  used  to  demonstrate  the  existence  of  general  positive  invariant 
sets . 

Proposition  4.1:  System  (2)  has  an  equilibrium  state  x^  4  0  if  and  only  if  for 

some  x  £  | (x  )  and  x  ^  0  either  Ax  =  B(x)  =  0  or  Ax  /  0,  B(x)  #  0  and  B(x)  £  I (A  x  ) . 
e 

Proof :  The  proof  follows  directly  from  the  fact  that  x^  is  an  equilibrium  state  if 

and  only  if  the  field  in  system  (2)  vanishes  at  x  .  Q 

e 

2 

If  on  the  boundary  of  an  open  set  in  R  ,  the  field  is  directed  everywhere 
toward  its  interior  then  that  set  is  positive- invariant.  Proposition  4.2,  an 
extension  of  the  well-known  Chetaev  Instability  Theorem,  establishes  this  fact 
formally. 

2 

Proposition  4.2:  Let  V  be  a  continuously  differentiable  functional  on  R  ,  let 
-1  + 

Q  4  <)>  separate  V  [  R  ]  and  let  x  =  f  (x)  admit  a  unique  solution  on  n.  If 
•-1  + 

ft  r-  V  [  R  ]  then  Q  is  a  positive  invariant  set. 


The  proof  follows  by  contradiction. 


\ 
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Proof :  Let  x^  e  ft  but  p(t;x^)  f  ft  for  some  t  t  IR  •  Since  V  is  continuous, 

-1  +  *  , 

V  [  R  ]  and  £2  are  open,  hence  t  =  inf{t  >  0|p(t;xo)  f  ft}  4  0  and 

*  -1  + 

p(t;x^)  e  ft  for  t  e  [0,t  ).  Since  ft  separates  V  [  R  ]  its  boundary  must  be  a 

-1  +  *  * 
boundary  of  V  [  R  ]  by  fact  A. 1  and  hence  V[p(t  Jxq)1  =  0.  Then  V[p(t  ;x  ) ] 

•k  • 

-  V[xQ]  =  -V(xQ]  <  0  and  for  some  t^  e  [0,t  )  we  must  have  V[p(t^;x  )]  <  0.  But 

•  —  2_  *4" 

this  contradicts  p(t^;x  )  c  ft  c  y  [  R  ].  □ 

b.  Some  Necessary  Conditions  for  Asymptotic  Stability  in  the  Large: 

From  the  results  of  Lyapunov  it  is  known  that  a  necessary  condition  for 
system  (2)  to  be  stable  is  that  the  spectrum  of  A  be  contained  in  the  closed  left 
half  of  the  complex  plane,  C~.  The  following  two  examples  indicate  that  system  (2) 
can  be  a.s.H.  even  when  A  is  a  stable  but  not  asymptotically  stable  matrix. 


Example:  4.1a: 


b: 


=  x2  +  .lx1(x1+x2) 

*2  '  ‘*1  +  -VvY 
h  ‘  'xi '  YYY 


(Figure  4a) 

(Figure  4b) 


x2 


X1(X1+X2) 


Hence  a  necessary  condition  for  a.s.£.  of  system  (2)  is  that  the  spectrum  of  A 
belongs  to  C-. 

If  the  stability  properties  of  the  linear  part  of  (2)  generally  determine  local 
behavior,  intuition  suggests  that  its  global  behavior  is  determined  by  the  quadratic 
part,  B(x).  This  is  partially  verified  by  the  following  three  lemmas. 

Lemma  4. 1:  If  | A |  f  0  and  B(x)  is  type  (i)  then  system  (2)  has  an  off-origin 
equilibr ium. 

Proof :  Since  x  t  0  implies  Ax  ^  0  and  B(x)  ^  0  under  the  hypothesis  above,  it 
suffices  to  show  that  for  some  x^  j*  0,  B(Xq)  e  |  (Ax^)  Proposition  4.1. 

Since  B(x)  e  | (Ax)  iff  |Ax,B(x) j  =  0, 


1 

f~ 

L  a2  J 


let  A  = 


(x.+x 


T 

a  x 
1 
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to  get  | Ax,B(x) I  = 


T 

a  x 
2 


T 

x  Gx 


T 

x  Hx 


=  x^qCx^/x^)  where  q  is  a  cubic  polynomial 


whose  leading  coefficient  is  al2h22  -  a22g22  =  a.  If  a  *  0,  there  exists  a 


real  root  v^  of  q(v)  =  0  and  hence 


some  B  e  R.  If  a  =  0  |AXq,B(Xq) |  = 


B 

-evo 
*12  S22 


a22  h22 


is  an  equilibrium  point  for 
0 


0  for  Xq  = 


1 


n 


Lemma  4.2:  If  | A |  =  0  and  B(x)  is  type  (i)  then  system  (2)  has  either  an  unstable 

origin  or  off-origin  equilibrium  states. 

T 

Proof :  Since  A  is  singular,  let  A  =  ab  .  If  there  exists  an  Xq  such  that 
0  4  Xq  i  l(b)  and  B(x^)  e  1(a)  then  by  Proposition  4.1  the  system  has  an  off-origin 
equilibrium  state.  If  for  all  x  4  0  B(x)  {  1(a)  choose  V(x)  =  a^x  to  get 
V (x)  =  a^[Ax  +  B(x) ]  =  a^B(x)  =  |B(x),a|  which  is  non-zero  for  all  x.  Since  B(x) 
is  continuous  it  follows  that  V(x)  is  sign  definite  and  since  V(x)  can  take  both 
positive  and  negative  values  in  the  neighborhood  of  the  origin,  the  latter  is  unstable. 

[In  [  2  ]  it  is  shown  that  when  B(x)  e  I  (a)  only  for  x  t  -t(b)  the  origin  is  unstable).  G 

T 

Lemmas  4.1  and  4.2  imply  that  B(x)  =  c  xDx  is  a  necessary  condition  for  global 
asymptotic  stability.  Following  the  results  of  section  3  one  might  expect  an 
additional  necessary  condition  requiring  D  to  be  focal.  The  following  two  examples 
demonstrate  that  system  (2)  can  be  a.s.St.  even  when  D  is  not  focal. 

Example  4.2: 

nr  lr  1  r- 

1  1 
0  1 


b. 


“  - 

[— 

— 

*1 

—i 

+1 

X2 

-i 

-1_ 

7*  “ 

- 

*i 

-i 

0 

X 

2 

0 

-1  _ 

+  x„ 


(Figure  5a) 


1 


'2  J 


+  x. 


0  0 


LX2J 


(Figure  5b) 


Figure  5a:  (Example  4.2a) 
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Examples  4.2a  and  b  are  a.s.i.  The  matrix  D  Is  critical  (i.e.  has  an  eigen¬ 
value  of  multiplicity  two)  in  both  cases  and  is  singular  in  case  b. 

When  D  »  db^  is  singular  B(x)  =  c^xdb^x  =  b^xdc^x.  In  such  a  case  there  is 
an  .tmbiguity  in  the  parametrization  of  B(x).  If  either  b  or  c  e  i(d)  then  the 
corresponding  B  matrix  is  critical  and  might  allow  a.s.Z.  behavior  as  in  example 
4.2b.  To  resolve  this  ambiguity  we  shall  classify  a  singular  matrix  D  as  nodal 
only  if  both  parametrizations  are  nodal  and  consider  only  the  critical  parametrization 
otherwise. 

T 

Lemma  4.3:  Let  B(x)  =  c  xDx  and  let  D  be  nodal.  Then  system  (2)  is  not  a.s.2. 

1’roof :  As  in  the  proof  of  Lemma  3.2  we  let  o(D)  =  { A , u }  X  i  |u|  with  no  loss  of 
generality  and  suppose  De  =  Xe.  The  lemma  is  proved  by  considering  the  two  cases 
Ae  e  | (e)  and  Ae  ^  | (e) . 

(i)  Case  Ae  e  1 (e) :  If  e  i  i(c)  then  Ae  ¥  0  implies  the  existence  of  an  off-origin 
equilibrium  by  Proposition  4.1.  If,  however,  Ae  =  0  then  since  e  is  an  eigenvector 
of  D  the  equation  (2)  has  a  finite  escape  trajectory  along  | (e)  as  shown  in  [1  ]. 

If  e  e  i(c)  then  Proposition  4.1  does  not  apply.  Using  the  change  of  basis 
[ , c ]  ^  we  obtain 

y 2  =  y2  (a22  +  d22y2^  where  d22  ^  °  ^since  D  is  noda1^  and  the 

system  is  unstable. 

T 

(ii)  Case  Ae  ^  ] (e) :  Choose  J  such  that  e^Ae  >  0. 

A  T 

From  fact  A. 4  choose  P  =  [e  (e  -ee)  ]  such  that 

Li  S 

c  i  C  (P)  and  c  (P)  C  C  ([PD]  ).* 

i  0  o  ~  o  s 


* 

The  assumption  that  D  is  nodal  is  crucial  here.  If  D  is  critical  and 
non-singular  then  Cp(P)  c  Cg([PD]g)  on^y  f°r  one  orientation  of  J  and  may  imply 

T 

Ae  <  0  (see  Appendix  A).  Hence  the  critical  case  is  considered  separately. 
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1  T  + 

Define  V(x)  =  —  x  Px  -  Y  where  Y  g  |R  and  Is  chosen  as  shown  below. 

V (x)  =  xTPAx  +  cTxxTPDx  and  for  x  e  C  ([PD]  ),V(x)  ^  0  if  and  only  if  cTx>  — ^ — 

°  S  xTPDx 

At  t  t 

where  a(x)  =  -x  PAx.  Since  a(e)  =  -£e  ee^Ae  <  0  (by  the  choice  of  P)  ®  assumes 


negative  values  on  -L(e)  and  by  Fact  A. 5  we  may  choose  Y  such  that  for  all 

x  e  H  =  C  (P)"H  (C)  -  C  [ PD]  ), 

Y  0  “  1)  s 

cTx  >  and  ii  c  V  IR+].  But  S2  also  separates  C  (p)  =  V  R+]  by  Fact  A. 3. 

x  PDx  "  Y 

Hence  fi  is  positive-invariant  according  to  Proposition  4.2,  and  since  0  k  the 

result  follows.  [1 

Since  Proposition  4.1  establishes  conditions  for  the  existence  of  off-origin 

T 

equilibria  depending  upon  the  zeros  of  |Ax,B(x)|  we  might  expect  B(x)  =  c  xDx  to 
necessitate,  in  turn,  0  4  |Ax,Dx|  =  |d||d  ^Ax,x|  (when  |d|  4  0)-  i.e.,  that  D  ^A  be 


focal.  The  following  are  examples  of  systems  with  D  A  not  focal  which  are  a.s.2. 


Example  4.3:  a. 


-1  1 


■1  0 


-1  0 


0  0 


0  -1 


1  0 


x  (Figure  6) 


x  +  (x1+x2) 


0  -1 


1  0 


x  (Figure  4b) 


The  matrix  D  A  is  critical  in  both  cases,  and  singular  in  4.3b. 

T 

Lemma  4.4:  If  B(x)  =  c  xDx  where  D  is  not  nodal  then  system  (2)  is  not  a.s.J.. 
when  both  |a|  =  | D |  =  0  or  either  A  ^D  or  8  ^A  are  nodal. 

2  -1  -1  i 

Proof:  e  e  R  is  an  eigenvector  of  D  A  or  A  D  iff  B(e)  £  I (Ae)  hence  if  either 

matrix  is  nodal  we  have  B(x)  e  | (Ax)  on  two  distinct  lines  | (e^)  4  | (e^) • 

Case  (i)  |d|  4  0:  Either  e.  or  e„  i  l(c) ,  hence  if  Ae.  4  0  and  e.  A  i(c)  then  the 

- -  1  l  l  j  T 

system  has  an  off-origin  equilibrium  according  to  Proposition  4.1. 

-1  T 

If  simultaneously  Ae^  =  0  and  e^  e  i(c),  that  is  D  A  =  c^b  for  some 


F igure 


\ 
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T  AT  •  t  T 

b  f.  i(e^))then  A  =  Dc^b  .  Defining  V(x)  =  x  [Dc^]^  implies  V(x)  =  c  xtDCj^lj^Dx 

T  T  T  T  i  i  T  2  •  ii 

=  (c  x)c  [J  D  JD]x  =  I D  |  ( c  x)  from  Fact  B.2.  Thus  V  ^  0  (since  |D|  i-  0  implies 

I D  |  >  0  for  non-nodal  D)  with  equality  only  for  x  e  i(c).  If  e ^  not  an  eigen¬ 
vector  of  A  then  a  complete  half  trajectory  cannot  remain  on  i(c),  hence  by 
Krasovskii's  extension  of  Lyapunov's  First  Instability  Theorem  [  3]  the  origin 
is  unstable.  If  Ae,,  t  |  ( e ,? )  then  D  =  “I  and  the  field  is  purely  quadratic  on 
l(e^)  in  the  direction  !(e^),  hence  the  origin  is  unstable. 

Case  (ii)  |p|  =  0:  Since  a  singular  non-nodal  matrix  is  x-critical,  we  have 

D  =  dd^  for  some  d  e  (R“.  If  Ad  i  j (d)  then  the  proof  is  identical  to  that  of 

case  (ii)  in  Lemma  4.3.  If  j  A |  =  0  and  Ad  e  | (d)  then  either  Ad  =  B(d)  =  0  or 

2 

Ax  t  | ( B ( x) )  for  all  x  e  R  with  off-origin  equilibria  in  either  case  according 

-1  —1  T 

to  Proposition  4.1.  If  | A |  ¥  0  and  Ad  c  | (d)  then  A  d  e  j (d)  and  A  D  =  ydd^ 

(y  f.  R)  .  Since  this  is  not  a  nodal  matrix,  it  need  not  be  considered  in  this 
proof.  LI 

In  direct  consequence  of  Lemmas  4.1-4  we  may  state  the  following  necessary 

conditions  for  a.s.H.  oc  (2). 

Proposition  4.3:  The  following  conditions 

(i)  A  has  eigenvalues  in  C- 
T 

(ii)  B(x)  =  c  xDx 

(iii)  D  is  not  nodal 

(iv)  either  A  or  D  A  exists  and  is  not  nodal 
are  necessary  for  system  (2)  to  be  a.s.ft. 

c .  Some  Sufficient  Conditions  for  Asymptotic  Stability  in  the  Large: 

According  to  Propositior  4.3  system  (2)  must  be  of  the  form 


x 


Ax  + 


T  ~ 
c  xDx 


(6) 


where  A  is  stable,  D  is  not  nodal  and  D  A  or  A  D  is  not  nodal.  In  this  section  we 
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derive  some  sufficient  conditions  for  a.s.i.  of  (6)  under  the  more  restrictive  assump¬ 
tion  that  D  is  focal. 


Corollary  4.1:  If  D  is  focal  then  (6)  has  no  off  origin  equilihrium  states  if  and 
only  if  D  ^A  is  either  focal  or  x-critical  where  x  r  l(c)  when  | A j  4  0  and  x  ^  i(c) 
when  j  A  j  =  0. 

Proof :  By  Lemma  4.4  D  ^A  must  be  either  focal  or  critical.  If  it  is  critical  but 

-1  i  2 

not  x-critical  then  D  A  =  «I  and  B(x)  e  ) (Ax)  for  x  c  R  with  resulting  off-origin 
equilibrium  states.  If  D  ^A  is  x-critical  then  B(y)  c  | (Ay)  if  and  only  if  y  t  | (x) . 
Under  this  condition  |a|  =  0  if  and  only  if  Ax  =  0.  Since  B(x)  =  0  if  and  only  if 
x  ■  i(c)  we  require  x  ^  i(c)  when  Ax  =  0  and  x  e  i(c)  when  Ax  ^  0  to  insure  against 
an  off  origin  equilibrium  state  according  to  Proposition  4.1.  f] 

Example  4.4  presents  a  system  which  satisfies  the  necessary  conditions  given 
in  Proposition  4.3  with  D  and  D  ^A  chosen  to  be  focal  matrices.  Since  this  system 
has  unbounded  solutions  it  is  apparent  that  further  restrictions  must  be  placed  upon 
A  and  D  to  obtain  sufficient  conditions  for  a.s.i. 

Example  4.4: 


“*  .  - 

— 

r~ 

r' 

x. 

-i  +i 

X, 

2 

-1 

x, 

1 

1 

+ 

1 

-i  -i 

X„ 

1 

1 

+2 

xo 

L.  2^ 

^  J 

L  2  - 

_ 

_ 

L  2J 

(Figure  7) 


(7) 


Along  the  line  x^  =  0  the  quadratic  part  of  (7)  vanishes.  Hence  in  the  vicinity  of 
this  line  the  influence  of  the  linear  part  is  felt  for  arbitrarily  large  values  of 
J|x)|  .  Considering  the  field  along  the  lines  x^  =  1  and  =  1  +  c  (e>0)  for  large 
positive  values  of  x^  we  have 

■  •  2 

Xj  =  1  and  Xj  =  (1+t)  -  ex^ 

x  2  =  x2  x2  =  (l+£)e  +  (1+2e)x2 


respectively.  Note  that  x^ 
Hence  along  the  strip  1  <  x^ 


increases 
<  1  +  e 


in  both  cases  while  x^  changes  direction. 


and  x„  >>  —  solutions  grow  in  an  unbounded 

2  c 


fashion. 


>  > 


1  is  crucial  to 


In  this  example  the  behavior  of  the  field  along  the  line  x^  = 
the  existence  of  unbounded  solutions.  Lemmas  4.5  and  4.6  and  their  consequences 
address  the  significance  of  this  affine  line  in  the  general  case. 

Lemma  4.5:  For  any  pair  of  matrices.  A,  D  £  and  any  x  t  R^,  if  a(x)  =  -Iax,x|, 

,  A  |  i  A  |  i 

°(x)  =  |Dx,x]  and  o(x)  =  -|Ax,Dx|  then  x  is  an  eigenvector  of  the  matrix  6(x)A  +  a(x)D 
with  corresponding  eigenvalue  n(x) . 

Proof :  Since  |(6(x)A  +  ot(x)D)x,  x|  =  6(x)|Ax,x|  +  a(x)|Dx,x|  =  0  we  have 

T  T 

[6(x)A  +  a(x)D]x  =  n(x)x  for  some  real  valued  function  n(x) .  But  x  xp(x)  =  x  [ 6 ( x ) A 
I  ^  T  T  T 

+  a(x)D]x  =  |  , tAx  Xt  =  [x,J  xj  [Ax , Dx]  =  -x  x|Ax,Dx|  hence  n(x)  =  -[Ax.Dx1. 

'x  JAx  x  JDx  _2 

Lemma  4.6:  The  quadratic  form  |Ax,Dx|  is  sign  definite  if  and  only  if  D  A  exists 

and  is  focal,  and  semi-definite  if  and  only  if  either  D  V  or  A  exists  and  is 


x-cr itical. 

Proof :  |ax,Dx|  is  sign-definite  if  and  only  if  it  is  never  zero  for  x  4  0,  in 

which  case  |  D  [  ^  0  and  |Ax,Dx|  =  |d|  |d  ^"Ax,x|  /  0,  which  is  equivalent  to  the 
statement  that  D  A  exists  and  is  focal.  Since  it  is  a  quadratic  form  |Ax,Dx[  is 
semi-definite  if  and  only  if  it  vanishes  on  exactly  one  line  through  the  origin. 

In  this  case  either  | A |  ^  0  or  |d|  4  0  (otherwise  the  form  would  vanish  on  at  least 
two  lines)  and  |Ax,Dx|  =  |d|  |d  ^Ax,x|  =  0  on  one  line  uniquely  is  equivalent  to  the 
statement  that  D  ^A  exists  and  is  x-critical,  with  a  similar  statement  for  | A |  #  0.  0 
Corollary  4.2:  If  D  is  focal  and  D  ^A  is  focal  (x-critical)  then  there  exists  a 


continuous  functional  u(x)  such  that  the  pencil  A  +  u(x)D  takes  eigenvalues  either 

-  +  T  ~+  2 

in  (R  or  iR  (  IR  or  IR  )  exclusively,  and  takes  x  as  an  eigenvector  for  all  x  e  R  . 


Proof :  If  D  is  focal  then  5(x)  =  |Dx,x|  ^  0,  hence  defining  u(x)  =  implies 

p(x)  is  continuous  and  x  is  an  eigenvector  of  A  +  p(x)D  with  eigenvalue  X(x)  = 

6(x) 

by  Lemma  4.5.  If  D  A  is  focal  (x-critical)  we  have  n(x)  sign  definite  (semi-definite) 


from  Lemma  4.6.  [3 


In  system  (6)  the  quadratic  part  vanishes  along  i(c)  and  the  influence  of  the 
linear  part  is  felt  in  the  neighborhood  of  this  line  even  for  large  values  of  x  as 


\ 
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41 

in  example  4.4.  If  we  define  =  p(c^)  then  by  Corollary  4.2  the  field  along  the 

T  A 

affine  line  c  x  =  |i^  is  given  as  f(x)  =  M  x  where  =  A  +  p  D  is  linear  and  takes 

U  n(c, V 


A  ' ' _  i'  * 

as  an  eigenvector  with  corresponding  eigenvalue  A^  =  ^  .  Lemma  4.7  demon¬ 

strates  that  the  sign  of  this  eigenvalue  is  a  crucial  indicator  of  a.s.f..  behavior 
of  system  (6). 

L emma  4.7:  If  D  is  focal  and  D  is  focal  or  x-critical  then  system  (6)  is  not 
a.s.7.  when  A  +  p(x)D  takes  eigenvalues  in  IR+. 

Proof :  If  A  +  p(x)D  has  positive  eigenvalues  then  p^  cannot  be  zero  since  A  has 

negative  eigenvalues.  Moreover,  o(M„)  =  {A  .61  c  IR+  with  M„c  =  A  r  from  corollary 

(r  00-  OiOi 

T  T  i  T 

4.1  hence  c  =  6^c  from  Fact  B.l.  On  the  affine  line  {x|c  x  =  p^}  we  have 

T  2 

f (x)  =  M^x  =  Mq(PqC  +  nc^)  n  e  E  and  hence  p^c  f(x)  =  p^fi^  2  0.  If  >  0  then  the 
half -plane  separated  from  the  origin  is  positive  invariant.  If  6^  =  0  then  the 
affine  line  is  itself  invariant.  In  either  case,  the  system  is  not  a.s.£.  □ 
Conversely,  it  is  easy  to  show  that  system  (6)  is  a.s.£..  when  the  pencil 
A  +  p(x)D  is  stable  if  we  insure  against  off-origin  equilibrium  states  according  to 
Corollary  4.1. 

Proposition  4.4:  The  following  conditions 

(i)  A  has  eigenvalues  in  C-I 
T 

(ii)  B (x)  =  c  xDx; 

(iii)  D  is  focal  and  D  A  is  either  focal  or  x-critical  where  x  e  i(c)  iff  | A |  4  0; 

(iv)  A  +  p(x)D  has  eigenvalues  in  IR  ; 
are  sufficient  for  (2)  to  be  a.s.5,. 

This  result  is  obtained  by  constructing  an  arbitrarily  large  neighborhood  of 
the  origin  and  demonstrating  that  it  is  positive-invariant.  Since  D  is  focal,  any 
trajectory  of  the  linear  system  (5)  -  x=Dx  -  contains  a  spiral, A,  corresponding  to 
one  period  of  oscillation  whose  endpoints  fall  on  the  same  side  of  the  same  line 


Here  and  in  the  sequel  we  assume  with  no  loss  of  generality  that  J|c||  =  1  in 
system  (6) . 


V 
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through  the  origin  of  the  phase  plane.  Joining  these  endpoints  with  the  line  segment 
A  yields  a  closed  curve  surrounding  the  origin  defining  the  neighborhood,  N ,  as 
depicted  in  Figure  8.  The  conditions  of  the  proposition  imply  that  the  linear 
field  is  always  directed  toward  the  interior  of  N  on  A,  while  the  entire  field  is 
directed  toward  the  interior  on  A  irrespective  of  the  orientation  of  its  linear  part. 
Proof :  Let  o(D)  =  {a  +  if?  ,a  -  ifSlwhere  we  assume  a  $  0  by  the  proper  choice  of  c 
(as  in  Lemma  3.2)  with  no  loss  of  generality.  Choose  an  orientation  for  J  such 

that  xTJ^Dx  >  0.  Define  A  =  {pet^c|t  c  [0,  }  with  endpoints  pc  and  pe^^'c  =  pyc 

8  * 

for  some  y  e  [0,1]  since  a  <  0  and  ■£—  is  the  natural  frequency  of  (5).  Further 
define  A  -  { C  c | Ce [ YP  »P ] }  to  close  the  boundary  of  W  as  depicted  in  Figure  8. 

The  scalar,  p,  will  be  chosen  below.  The 

normal  vector  to  the  curve  A  at  any  point 

A 

x  is  given  as  nA(x)  =  JDx  and  this  normal 

is  directed  into  the  interior  of  N  since 

x^n^fx)  =  x^JDx  =  -xTJTDx  <  0.  Similarly, 

the  interior  directed  normal  on  A  is  given 

by  n^(x)  =  Jc.  It  now  suffices  to  show  that 
T  T 

n^(x)f(x)  and  n^(x)f(x)  are  non-negative  for 
x  e  A,  A  respectively.  The  first  inequality 

follows  since  x  c  A  implies  x  =  tc 

T  TT  2  TT  ±  ± 

r,  (.  [yp»p]  hence  n^(x)f(x)  =  c  J  [CAc  +  C  Dc]  =  t[c  J  Ac  +  CciJ'lDc]  which  is  always 

T  T 

positive  for  large  p  since  c  J  Dc  >  0.  To  see  the  second  inequality  holds  note 
n  (x)^f(x)  =  x^D^J^fAx  +  c^xDxJ  =  x^D^J^Ax.  But  x^D^j^Ax  must  have  the  same  sign  as 


Figure  8 
T.T 


T  T 

x  J  Dx  since 


T  T  T 
x  D  J  Ax 


T  T 
x  J  Dx 


Ax ,Dxj  _  j  q  from  Corollary  A. 2. 


I Dx,x I 

Since  p  may  be  arbitrarily  large  in  the  above  argument,  we  have  shown  that  all 

T  T 

solutions  to  (6)  are  bounded.  If  c  J  Ac  >  0  then  the  argument  holds  for  all  values  1 

T  T 

of  p  and  the  proof  is  complete.  If  c  J  Ac  <  0  then  it  is  easy  to  show  that  no  limit 

cycles  exist.  Since  the  conditions  of  this  proposition  exclude  off-origin  equilibrium 

states  as  well,  the  boundedness  of  solutions  implies  that  the  origin  is  a.s.fc.  □  y 
*  Note:  if  y  =  1  then  a  =  0  and  D  ijj,  skew-symmetric  in  some  basis  for  which  V=x  x 
Is  a  Lyapunov  function  -  i. e.  V  >  0,  V  =  x  Ax  <  0  under  these  hypotheses. 
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d .  Special  Cases  -  Critical  Matrix  D: 

Propositions  4.3  and  4.4  give  necessary  and  sufficient  conditions  respectively 

T 

for  a.s.£.  of  system  (2).  In  both  cases  B(x)  =  c  xDx  so  that  system  (2)  can  he 


described  by  (6).  The  necessary  conditions  stipulate  that  the  matrix  D  be  non- 


nodal  while  the  sufficient  conditions  in  Proposition  4.4  are  given  for  the  case 
when  D  is  a  focal  matrix.  To  complete  the  investigation  of  the  stability  properties 


of  (2)  we  consider  in  this  section  system  (6)  when  D  is  critical.  This  class  is 
most  easily  treated  by  dist inguishing  between  the  cases  when  D  is  non-singular  and 


D  is  singular. 

The  case  when  j  D |  4-  0  closely  resembles  that  analyzed  earlier  in  section  4c 
where  D  is  focal  but  D  *A  is  critical  in  that  the  stability  properties  of  both 


systems  are  seen  to  depend  on  the  vector  c. 

D  Critical  and  Non-Singular: 

Lemma  4.8:  If  D  is  critical  and  non-singular  then  a  necessary  condition  for  (6) 
to  be  a.s.H.  is  that  D  be  c  -critical. 

i 

This  Lemma  is  proved  by  demonstrating  the  existence  of  a  positive  invariant 
region  disconnected  from  the  origin.  If  V(x)  =  (x-m) "V (x-m)  -  Y,  for  some  constant 

vector  m,  matrix  P,  and  scalar  y ,  then  a  separation  of  the  region  V(x)  >  0  can  be 


shown  to  be  contained  in  the  region  V(x) 
a.s.I.  according  to  Proposition  4.2. 
Proof :  If  D  is  critical  but  not 

2 

c  -critical  there  exists  a  vector  e  e  R 

L 

T 

such  that  De  =  Xe,  X  >  0  and  c  e  >  0. 

j  T 

From  Fact  A. 4  choose  P  =  [e,(e  -£e)  ] 

i.  J.  3 

such  that  c^  |  C^(P)  ‘  ^qOPD]^)  and  we 

T 

may  choose  m  i  Cq(P)  with  c  m  >  0.  De¬ 
fine  V(x)  =  l/2(x-m)TP(x-m)  -  y  for 
some  Y  c  R  to  be  chosen  below.  The 


0  demonstrating  that  the  origin  is  not 


C0([PD],)+m 


Figure  9 


■v 
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region  V(x)  >  0  is  depicted  in  Figure  9  and  contains  the  open  set  ft  consisting  of 

-1  +  T  -1  + 

those  elements  x  £  V  [  R  ]  for  which  c  x  >  0.  Formally,  C  (p)  +  m  =  V  [  |R  ] 

and  Si  =  (C  (P)nH^(c)]+  m  separates  V  IR+]  as  shown  in  Fact  A. 3.  If  y  =  x  -  m 

then  x  £  V  \  R  ]  iff  y  £  C  (P).  Since  V(x)  =  (x-m)^Px  =  y^Pv  we  have  from  (6) 

V(y)  =  cTy  yTPDy  -  a(y) 

T 

where  a(y)  is  a  functional  with  linear  and  quadratic  terms.  Since  y  PDy  >  0  when 

v  t  C  (P) ,  V(v)  >  0  on  that  set  iff  c^y  >  ° .  Since  a(y)  can  be  shown  to  be 
Y  '  y  PDy  T 

negative  on  the  common  boundary  of  C  (P)  and  C  ([PD]  )  in  the  half-plane  c  x  >  0  - 

0  Os 

i.e.  on  ((e) nH^ ( c )  -  it  follows  from  Fact  A. 5  that  a  Y  e  R  exists  such  that  the 
inequality  holds  for  y  e  C^(P)nH^(c).  This  is  equivalent  to  the  statement  V(x)  >  0 
for  x  e  ft. 

Since  ft  separates  V  p  ]  and  V(x)  >  0  in  0,  SI  is  positive-invariant  by 

Proposition  4.2.  Since  0  ^  ft  the  system  is  not  a.s.i.  0 

Lemma  4.9:  If  D  is  c^-critical  and  non-singular  a  necessary  condition  for  (6)  to 
be  a.s.^.  is  that  D  A  be 
either  (i)  focal 

or  (ii)  singular  and  x-critical  where  x  4  J-(c). 

Proof :  By  Lemma  4.4,  (6)  is  not  a.s.i.  when  D  A  is  nodal.  The  result  follows  by 

showing  that  the  cases  where  D  ^A  is  critical  and  either  non-singular  or  c^-critical 
are  not  a.  s.  #■ . 

If  D  ^A  is  critical  and  non-singular  with  eigenvector  x  ^  i(c)  then  B(x)  e  |(Ax), 

B(x)  4  0  and  Ax  4  0.  Hence  an  off-origin  equilibrium  exists  by  Proposition  4.1. 

If  D  ^A  is  c^-critical  then  c^  is  also  an  eigenvector  of  A  since  D  is  c^-critical. 

A  -1 

Under  the  change  of  basis  y  =  [c^c]  x  equation  (6)  yields  y^  =  y2^22^2+322^ '  Since 
^22  ^  ®  the  system  has  unbounded  solutions.  G 

Corollary  4. 3:  Let  D  be  c^-critical  with  eigenvalue  Y  <  0.  If  D  \  is  focal  or 

x-critical  then  the  pencil  A  +  p(x)D  has  eigenvalues  in  IR  for  all  x  {  i(c)  iff 
T  T 

r  Ac  x  D*  1  0. 

i  l 
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-  |  Ax  f  Dx  | 

Proof :  For  x  ^  ±(c)tA  +  p(x)D  is  defined  and  has  an  eigenvalue  X(x)  =  -r - r— 

|  Dx , x | 

according  to  Lemma  4.5  where  X  is  sign  definite  (or  semi-definite)  according  to 

T 

Lemma  4.6.  If  |Dx,x|  =  x^Dx  >  0  then  X(x)  <  0  iff  |Ax,Dxj  5  0.  Since  |Ac^,Dc^| 

1  1  T 

=  y|Ac^,c^|  =  -yc  Ac^  and  y  <  0  the  result  follows,  f] 

These  results  lead  to  the  following  necessary  and  sufficient  conditions  for 
a.s.J.  when  D  is  critical  and  non-singular. 

Proposition  4.5:  If  D  is  critical  and  non-singular  then  system  (6)  is  a.s.Z.  if 

and  only  if 

(i)  D  is  c  -critical 

(ii)  D  ^ A  is  either  focal  or 

D  ^A  is  singular  and  x-critical  where  x  ^  i(c) 

2 

(iii)  For  some  x  c  H  the  pencil  A  +  p(x)D  has  a  negative  real  eigenvalue. 

Proof :  As  in  the  proof  of  Proposition  4.4  we  may  assume  that  Dc  =  yc  where  y  <  0 

T 

and  x^Dx  >  0  with  no  loss  of  generality.  Hence,  under  the  conditions  (i)  and  (ii) 

T 

condition  (iii)  is  equivalent  to  the  condition  c  Ac  >  0  by  Corollary  4.3,  or 
T 

''DxJ^Ax  >  0  by  Lemma  4.5  and  4.6. 

Necessity:  By  Lemma  4.8  and  4.9  we  require  (i)  and  (ii).  If  (iii)  does  not  hold 
T 

then  c  Ac  <  0  as  shown  above.  In  this  case  the  proof  is  identical  to  that  in 
Case  (ii)  of  Lemma  4.3,  where  (6)  is  shown  to  be  not  a.s.H.  by  the  construction  of 
a  positive-invariant  set,  C  (P) ,  disconnected  from  the  origin. 

Suf f iciency:  is  shown  by  constructing  an  arbitrarily  large  Jordan  curve  f  as  depicted 
in  Figure  10, such  that  the  field  at  every  point  on  the  curve  is  directed  towards 
the  interior  of  the  enclosed  region  N.  Since  D  is  c^-critical  the  curve  etDx, 
t  t  R,  x  e  H^(c)  intersects  |  (y)  for  all  y  ^  i(c)  (Figure  10a).  The  curve  y  may 
then  be  expressed  as  the  union  (Figure  10b) 

-  —  +  + 

f  5  A  uA  uA  ua  ur 
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where  A  and  A  are  segments  of  curves  of  the  form  e^x  and  A  ,A+  and  T  are  line 
segments,  all  of  which  are  defined  as  follows. 

A  =  {n(x0~tb)|t  e  [O.t^]};  A  =  {pe  tDx^|t  e  f 0 , t ^ j 

A+  =  (n(x2+tb)|t  e  [0,t3]);  A+  =  (ne  tDx3|t  e  [0,t^]l 
and  r  =  { nfx^-yc)  |u  e  [O.y^-Ug]} 


where  xQ  A  -Cj>  +  V  ;  b  ±  c±  +  £qC 


£0>  ° 


A 

\in  =  max 


c^Ac  cJac± 


T  *  T 

ciDc  CiDci 


n  >>  1  and 


X1,X2,X3  an<*  x4  are  enc*  P°^nts  °f  A  ,A  ,A  +  and  A+  respectively.  The  parameters 

t^  and  t^  are  chosen  so  that  A  and  A  cross  i(c),  and  t^  is  chosen  so  that  e  |  (c^ 

-  c,c)  where  0  <  c  <  e_.  The  parameters  t  and  u,  are  such  that  the  curve  is  closed 
1  10  A  4 

at  x,  . 


A  family  of  such  Jordan  curves  can  be  generated  by  changing  the  parameter  n  and 

the  arguments  given  below  apply  to  all  members  of  the  family  for  n  sufficiently  large. 

As  in  the  proof  of  Proposition  4.4  the  field  is  directed  into  the  interior  of  N 
+  T 

on  A  and  A  since  (Dx)  Ax  >,  0  as  shown  above.  The  interior  normal  direction  on  T 
T  T 

is  c  and  c  f(x)  >  0  since  c  Dc  >  0  and  the  quadratic  term  dominates  the  linear 
it  l 

term  of  the  field  by  the  choice  of  On  A  and  A  the  linear  field  dominates,  and 

T 

is  directed  towards  the  interior  of  N  since  c  Ac^  >  0  as  shown  above,  q 
D  Critical  and  Singular: 

Proposition  4.6:  If  D  is  critical  and  singular  then  system  (6)  is  a.s.I.  iff  | A |  ?  0 
and  A  =  yD  for  some  y  e  R,  in  which  case  the  pencil  A  +  y(x)D  has  eigenvalues  in 
R  when  it  is  defined. 

Proof :  From  Lemma  4.4,  if  If)}  =  0  then  |  A  |  0  and  A  D  is  not  nodal.  Since 

T  2-1  —  j  -  -] 

D  =  dd^  for  some  d  e  |R  ,  A  D  is  nodal  unless  A  d  r  | (d)  in  which  case  A  D 


YddJL  *  yD  for  some  Y  £  R. 


-1 


-1. 


1 

"an 

1  ■ 
a12 

'  0 

1~ 

we  have  A  * 

0 

a22 

J 

,  D  » 

0 

0 

hence  y *  a22Y2  and 


Figure  10 
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a22t 


a22fc 


a22t 


y2(t)  =  e  y20  "  °*  In  thls  case  yi  =  (ail+Cie  y20)yl  +  a12e  y20  +  c2e 
and  y^(t)  0.  The  system  is  a.s.H. 

In  the  y  coordinate  system,  A  +  u(y)D  has  the  eigenvalue 

2 


2a 22 1  2 


20 


A  (y) 


.  -  lAy.Dyl  _  a22y2  ,  Q 

|fty.y|  yl 


according  to  Lemma  4.5.  0 

e.  Summary  of  Results: 

Sections  4  a-d  have  investigated  the  properties  and  conditions  for  a.s.H.  of 
system  (6)  according  to  the  scheme  depicted  in  Figure  6.  By  Proposition  4.3  in 
section  4b  matrices  D,  A  D  or  D  A  cannot  be  nodal  if  system  (6)  is  to  be  a.s.H. 


A  stable,  B(x)  =  c  xDx 


In  section  4c,  Proposition  4.4  specifies  sufficient  conditions  for  a.s.i.  which  are 
also  seen  to  be  necessary  when  applied  strictly  to  the  right  branch  of  the  tree  in 
Figure  11, where  D  is  focal.  In  section  4d  Propositions  4.5  and  4.6  specify  sufficient 
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conditions  for  a.s.£.  which  are  also  necessary  when  applied  to  the  left  branch  of 
the  tree  where  D  is  critical.  The  stability  condition  on  the  pencil  A  +  M(x)D  is 
seen  to  have  the  intuitively  appealing  geometric  interpretation  that  the  perturbation 
of  the  pure  quadratic  field  by  the  addition  of  the  linear  field  must  be  directed 
toward  the  origin  at  every  point.  The  results  of  sections  4  a-d  lead  to  the 
central  contribution  of  this  paper  given  by  Theorem  2. 

Theorem  2:  The  following  conditions  are  necessary  and  sufficient  for  a.s.i.  of 
system  (2) : 


(i)  the  eigenvalues  of  A  lie  in  C 


(ii)  B(x)  =  c  xDx 

(iii)  D  and  D  \  or  A  exist,  are  either  focal  or  critical,  and  satisfy 
one  of  the  following: 

(a)  D  and  D  ^A  are  focal; 

(b)  D  is  focal  and  D  *A  is  x-critical  where  x  e  i(c)  iff  | A |  i  0; 

(c)  D  is  c^-critical  and  non-singular,  and  D  ^A  is  focal; 

(d)  D  is  c^-critical  and  non-singular,  D  ^A  is  singular  and  not  c^-critical; 

(e)  D  is  singular  and  A  =  yD  for  some  scalar  y. 


(iv)  if  y (x)  = 


Ax,x 


Dx,x 


then  the  pencil  A  +  y(x)D  has  a  negative  real  eigenvalue 


for  some  x  e  IR  . 

Proof :  Conditions  (i),  (ii),  and  the  general  statement  of  (iii)  are  necessary 
according  to  Proposition  4.3.  Under  these  conditions  (iv)  is  equivalent  to  the 

2 

statement  that  A  +  y(x)D  has  negative  real  eigenvalues  for  all  x  e  IR  where  y  exists 
according  to  Lemmas  4.5  and  4.6.  If  D  is  focal  then  the  additional  condition  (iiia) 
or  (iiib)  is  necessary  and  sufficient  according  to  Proposition  4.4.  If  D  is  non¬ 


singular  and  critical  the  additional  condition  (iiic)  or  (iiid)  is  necessary  and 
sufficient  according  to  Proposition  4.5.  If  D  is  focal  then  ( i) , ( ii) , ( iiie)  are 
necessary  and  sufficient  according  to  Proposition  4.6.  □ 
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Two  examples  of  systems  satisfying  the  conditions  of  Theorem  2  are  given 


below  with  accompanying  figures. 


Example  4.5: 


—  — 

r  — 1 

-1  0 

0  -1 

a.  x  = 

X  +  X, 

0  -1 

1 

1  0 

_  J 

_  _ 

~-l  2 

r 1  n 

b.  x  = 

X  +  X 

-2  -1 

1 

-i  i 

— 

_ 

(Figure  12a) 


(Figure  12b) 
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Appendix  A 

2 

Some  Geometrical  Results  In  IR 

Many  of  the  apparently  complicated  proofs  in  the  body  of  this  paper  are  based 
upon  some  simple  geometrical  notions.  Unfortunately  the  algebra  and  notation  needed 
to  make  these  concepts  precise  are  unwieldy.  This  appendix  presents  the  most  im¬ 
portant  of  these  ideas  as  Facts  A. 1-5  of  which  A. 4  and  A. 5  are  proved  in  detail 

since  they  are  used  throughout  the  paper. 

2 

Def inition  A. 1:  If  9  c  IR  is  open,  the  boundary  of  6  is  defined  as 

A  - 

36  =  8  -  e 

2 

Fact  A. 1 :  Let  6  c  r  be  open  and  separated  Ly  6^.  Then  30^  c  30.  If  =  0 

and  0^002  =  4>  this  merely  states  that  =  <t>- 

Def inition  A. 2:  Let  e,e  be  linearly  independent  in  R  .  Then 

C  (e,e)  =  {ae  +  6e|a0  >  y  >  0} 

Y 

CQ(e,e)  defines  a  conical  region  in  R2.  If  >  y^  C ^  (e,e)  is  a  proper  subset  of 
C  (e,e) .  The  distance  of  the  set  C  (e,e)  from  the  origin  increases  monotonically 

y2  Y 

with  y  and  tends  to  infinity  as  y  tends  to  infinity. 

An  alternate  description  of  CQ(e,e)  and  Cy(e,e)  uses  and  e^. 


x  c  C_(e,e)  iff  -x  e  eTx  >  0 

0  j  i 

'K  ,  r  vT 


x  c  C  (e,e)  iff  -x  e  e  x 
Y  11 


T'l  2 
feiei] 


>  Y 


If  P  =  -[eel  then  C  (e,e)  =  C  (P)  defined  in  section  3a. 

1  1  s  Y  y 

7  T'v.  .1 

[e  e  ] 

1  i 

'v.  2  %  . 

Fact  A. 2:  Let  e,e  and  x  e  IR  .  If  e  and  x  lie  in  the  same  half  plane  defined  by  | (e) 


I  r  ^ 

and  x  $  Cg(e,e)  then 


cQ(e,e)  -  c0(e,x) ' 
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Fact  A. 3:  Let  c  A  C  (e,e).  If  H  (c)  =  {x|cTx  >  n)as  defined  in  section  3a,  given 
-  x  0  n 

an  n  Y  0  there  exists  a  y  >  0  such  that 


H  (c)nC  (e,e)  separates  C  (e,e) 

n  y  y 

If  c^  e  3C^(e,e)  then  the  above  result  holds  only  for  n  =  0. 

Fact  A. 4:  Let  D  be  critical  or  nodal  with  o(D)  =  {X,y}  and  X  >  |y|.  Then  given 

any  c  e  B  there  exists  an  indefinite  matrix  P  such  that  c^A  Cq(P)  ant*  Cq(p)  - 

Comment:  Fact  A. A  is  used  throughout  the  paper  to  prove  that  the  given  system  is  not 

2 

a .  s .  £ .  The  principal  idea  is  to  establish  a  region  in  (R  by  the  proper  choice  of 
P,  which  is  positive  invariant.  Embedding  C  (P)  in  C  ([PD]  )  establishes  the 
positive  invariance  of  one  half  of  the  cone  C^(p)  under  the  influence  of  the  pure 

quadratic  field.  Fact  A. 4  states  that  if  D  is  nodal  or  critical  such  an  embedding 

is  always  possible. 

^  A  n. 

Proof :  Let  De  =  Xe  and  define  e  =  e  +  £e^  for  some  £  >  0.  Note  ex  =  ex  “  Ee* 

A  „  VT. 

Define  P  =  -[e  e  J  to  get 
1  i  s 

™  1.  TtIl  1,  <vT_  'v  T. 

PD  =  - [e  e  +  e  e  ]D  =  -  — [e  e  D  +  ye  e  I 

2  x  x  i  x  2  x  i  Xi 

T  T 

since  e  D  =  ye  from  Fact  B.l.  Thus 
x  x 

[PD]  =  -[e  gT] 
s  is 

where  g  =  -|[DT  +  yl]e^.  If  g  e  x(e),  then  it  directly  follows  that  0  =  eTg  =  -c/2 
T 

(X+y)e  e,  hence  X  =  -y  and 


CQ([PD]s)  -  R  -  |(e)  a  CQ(P). 


If  g  i  x(c)  we  consider  the  three  cases  where  D  is  (i)  critical  but  not  x-critical 
(ii)  x-critical  and  (iii)  nodal. 

In  case  (i)  D  =  oI(a  >  0)  and  PD  =  aP  so  that  CQ(P)  =  CQ([PD]s).  For  cases 
(ii)  and  (iii)  we  demonstrate  that  CQ(P)  c  CQ ( [PD] g)  by  the  equivalent  demonstration 


that 
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Ve*e>  1  Co(e,_8i) 


using  the  results  of  A. 2. 

T  T  T  ^ 

Since  eg  =  e  g  =  -e/2(A+u)  <  0  and  e  e  =  ee  e  >  0  it  follows  that  e  and 
-L  1  111 

-g^  lie  in  the  same  half  plane  defined  by  | (e) .  To  show  -g^  |  CQ(e,e)  it  suffices 
T  T  'vT  T 

to  show  0  >  g^Pg^  =  “(g^^)  (eigi)  •  Since  g^e^  <  0  from  above,  this  follows  if  and 

vr  t  eiDe  x 

only  if  e  g  =  -e/2  e  e(A-y  +  e  — ^ —  )  <  0.  The  last  inequality  is  evidently  true 

e  e 

for  small  enough  e  if  D  is  nodal  since  A  >  |y|.  If  D  is  x-critical  then  A  -  y  =  0 

VT  A  * 

and  the  sign  of  e^g^  depends  upon  J  -  i.e.  the  orientation  of  e  =  Je.  In  either 

^T  i  „  % 

case  e1g1  <  0  and  g±  $  C  (e,e).  Thus,  by  Fact  A. 2 


CQ(P)  =  CQ(e,e)  =  C0(e,-8i)  =  C0([PD]s). 


i  t  t  2  T  T 

To  show  that  cj:  CQ(P)  we  note  that  c^Pc^*  -(c*e^)  +  tCc^e^e  O  f  0  for  small 

enough  e.  □ 

v  T  T 

Fact  A. 5:  Let  e,e,c,D  and  P  be  as  in  A. 4.  Let  a(x)  =  x  Qx  +  g  x  for  arbitrary 

T 

matrix  Q  and  vector  g  (possibly  0)  such  that  e  Qe  >  0.  Then  for  some  y  t  R,  there 

T  T 

exists  a  separation,  11,  of  C^(P)  such  that  T(x)  =  a(x)  +  c  x  x  PDx  >  0  for  x  £ 
Comment :  This  fact  is  used  to  prove  the  existence  of  positive-invariant  sets  dis- 
connected  from  the  origin.  If  x  =  b+Ax+c  xDx  as  in  section  4  and 
V(x)  =  -j  xTPx  -  Y,  then  V(x)  =  H^x). 

Proof :  If  x  e  Cq([PD]^)  then  'f(x)  >  0  iff  c  x  >  0(x)  where  9(x)  =  .  We 

&  x  PDx 

first  show  for  some  Y,sup  0(x)  =  6  <  «.  This  follows  since  j 0 (x ) |  <  <*  as 

xcC  (P) 

C^(P)  c  Cq([PD]s),  while  on  the  common  boundary 

8C0(P)n3C0({PD]s)  =  *U)’ 


choice  of 


T  T  T  A 

Unless  |  D j  =  0  in  which  case  D  =  ee^  and  e^D  e  =  )]  e]j  regardless  of  the 
if  J . 
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whtch  is  arbitrarily  close  to  3 C  (P) ,  we  have  a(x)  >  0  for  x  =  Be  when  B2  >  6n, 

Y  u  * 

T 

for  some  ^  e  R  since  e  Qe  >  0. 

Case  (i)  e  t  i(c):  By  Fact  A. 2  there  exists  a  Y  such  that  C  (P)nHrt  (c)  separates 

y  6() 


C  (P).  Hence  defining  ft  =  C  (P)ntfg  (c)  yields  the  desired  result. 

0  2 

Case  (ii)  e  e  i  (c) :  Since  a(x)  >  0  for  x  =  Be,  B  >  Bq  by  continuity  a(x)  >  0 
2 

for  x  =  Be'  8  >  8^  for  some  e'  in  a  neighborhood  of  e  contained  in  C  (p) .  Hence  if 


6  =  max( 80,B1,eo/cTe' ,60/c^e' } 


and 


%  5  >  «0  1  62  S  'T'’S  >  % 


AT  T 

then  a(x)  >  0  in  the  rectangle  n  =  (x|c  x£[O,02]  and  c  x  >  0^}.  Choosing  Y  such 
that  ft  =  C  (p)nHQ  (cx)  separates  CQ(p)  according  to  Fact  A. 2  we  have 


while 


T 

c  x  > 


0  > 


~a(x) 

T 

x  PDx 


for  x  £  ftnH 


cTx  >0  >  0  >  for  x  e  ft  -  n  as  shown  in  Figure  A.l.  □ 

1  0  x  PDx 

Comment:  When  i(c)  does  not  coincide  with  the  common  boundary  of  C  (P)  and  C  ([PD]  ) 
-  0  Os 

it  is  possible  to  separate  C  (P)  using  the  half  plane  H  (c)  as  shown  in  case  (i). 

'  '  Y  0Q 

When  i(c)  coincides  with  this  boundary,  the  proof  is  more  involved  as  shown  in 
case  (ii).  The  region  ft  is  expressed  as  the  union 


ft  =  (linn)  u(Hq  (c±)nft) 


-o(x) 

T 

x  PDx 


or  H’(x)  >  0. 


Figure  A.l 
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Appendix  B 

Some  Properties  of  (2x2)  Matrices 


Fact  B.l:  If  x  is  an  eigenvector  of  the  matrix  A,  x^  is  an  eigenvector  of  the 

T  T 

matrix  A  .  If  0(A)  =  {X,y}  c  IR  and  Ax  =  Xx,  then  A  x^  =  yx  . 

T 

Proof :  (Ax.x^)  =  (Xx,x±)  =  (x,A  x±)  =  0  where  (x,y)  is  the  inner  product  of  the 

T 

vectors  x  and  y.  Hence  x^  is  an  eigenvector  of  A  • 

T 

Let  the  eigenvalue  of  A  corresponding  to  x^  be  a  and  let  the  eigenvector  of 

T 

A  corresponding  to  the  eigenvalue  M  be  y.  A  x^  =  ax^ ;  Ay  =  uy. 

T 

(Ay ,x^)  =  (uy,x^)  =  ( y , A  x^)  =  (y.ax^) 

Hence 

T 

(y-a)y  x^  =  0 


Since  the  eigenvectors  y  and  x  are  independent  U  =?.  0 

!  0  -1 

Fact  B.2:  If  A  is  a  non-singular  matrix  and  J  = 

L1  ° 


then  A 


0  -1 

all  a12 

T.T, 

a22  “ai2 

,  A  = 

J  A  J  = 

1  0 

_a2l  a22  . 

_'a21  all_ 

1 

I A I 


JTATJ. 


.  □ 


□ 


Proof:  If  J  = 
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